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ABSTRACT 


Ferret systemic coronavirus (FRSCV) infection is associated with an emerging, highly fatal disease of 
ferrets. Enhanced macrophage tropism and the resulting induction of pyogranulomatous lesions are 
shared with feline infectious peritonitis virus (FIPV) infection in cats, but are not features of ferret enteric 
coronavirus (FRECV) infection. Comparative sequence analysis of the distal one-third of the genomes 
of one FRSCV and one FRECV strain showed that these two ferret coronaviruses share >96% nucleotide 
sequence identities in the membrane (M), nucleocapsid (N) and non-structural protein genes (partial 
polymerase, open reading frames [ORFs] 3 and 7b). The envelope (E) protein gene showed a moderate 
nucleotide sequence similarity of 91.6%. In contrast, nucleotide and amino acid sequence similarities 
observed with the spike (S) protein were only 79.5 and 79.6%, respectively. Twenty-one amino acid 
differences within a 195-199-amino acid C-terminal portion of the S protein were conserved between 
3 strains each of FRSCV and FRECV. Both systemic and enteric strains were found to carry a single ORF 
3 gene with truncated proteins observed in two out of three FRSCV strains examined. The two enteric 
strains analyzed each contained an intact ORF 3 gene. Phylogenetically, FRSCV is more closely related to 


FRECV than to other group 1 coronaviruses. 


© 2010 Elsevier B.V. All rights reserved. 


1. Introduction 


An emerging systemic disease in ferrets characterized by feline 
infectious peritonitis (FIP)-like clinical signs and lesions, has been 
recently observed in ferrets across the U.S. and also in Europe 
(Garner et al., 2008; Juan-Sallés et al., 2006; Martinez et al., 2006, 
2008). Clinical findings commonly observed in the FIP-like disease 
in ferrets include anorexia, weight loss, diarrhea and the pres- 
ence of large palpable intra-abdominal masses. Strikingly similar 
to FIP is the presence of widespread granulomas on serosal sur- 
faces and within the parenchyma of the abdominal and thoracic 
organs. Histopathologic examination shows a systemic pyogranu- 
lomatous inflammation involving the mesenteric adipose tissue, 
lymph nodes, liver, lungs, kidney, spleen, pancreas and adrenal 
glands. Moreover, the immunohistochemical detection of coron- 
avirus antigen in the cytoplasm of intra-lesional macrophages is 


* Corresponding author at: Diagnostic Center for Population and Animal Health, 
Rm. 136E, 4125 Beaumont Road, Lansing, MI 48909, USA. Tel.: +1 517 432 5794; 
fax: +1 517 353 5096. 

E-mail addresses: wise@dcpah.msu.edu (A.G. Wise), kiupel@dcpah.msu.edu (M. 
Kiupel), zoopath@aol.com (M.M. Garner), maes@dcpah.msu.edu (R.K. Maes). 


0168-1702/$ - see front matter © 2010 Elsevier B.V. All rights reserved. 
doi:10.1016/j.virusres.2009.12.011 


consistent with what has been observed in tissues of cats with 
FIP. 

FIP is a fatal, multi-systemic, immune-mediated disease of cats 
caused by FIP viruses (FIPV), postulated to be virulent mutants 
arising spontaneously from persisting low pathogenic or non- 
pathogenic feline enteric coronavirus strains (De Groot-Mijnes et 
al., 2005; Poland et al., 1996; Vennema et al., 1998; Herrewegh et 
al., 1995; Rottier et al., 2005; Chang et al., 2009; Pedersen, 2009). 
More recently, a study by Brown et al. (2009) called into ques- 
tion this in vivo mutation hypothesis of FIP pathogenesis. Their 
findings supported another hypothesis, namely, that genetically 
distinctive avirulent and virulent forms of feline coronavirus (FCoV) 
co-circulate in natural populations. 

In a recent publication (Garner et al., 2008), we reported the 
detection of coronaviral RNA in the tissues of the affected ferrets 
by amplification of a portion of the coronavirus spike gene using 
consensus RT-PCR. A unique 599-bp partial spike sequence of the 
ferret systemic coronavirus (FRSCV) was derived and analyzed. This 
sequence was found to be more similar (77% deduced amino acid 
sequence similarity) to the previously described (Wise et al., 2006) 
ferret enteric coronavirus (FRECV) than to other group 1 coro- 
naviruses, feline coronavirus, transmissible gastroenteritis virus 
(TGEV) of swine and canine coronavirus (CCV) (71-73% deduced 
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amino acid sequence similarities ). FRECV is associated with a differ- 
ent disease in ferrets, epizootic catarrhal enteritis (ECE) (Williams 
et al., 2000). ECE is a highly contagious diarrheal disease of ferrets 
associated with a mortality rate of <5%. When first recognized, the 
disease was more severe in older ferrets but often mild and subclin- 
ical in juveniles. We previously characterized FRECV, and found it 
to be a novel group 1 coronavirus (Wise et al., 2006), based upon 
analyses of partial sequences of the polymerase (pol), spike (S) and 
membrane (M) protein genes and the complete sequence of the 
nucleocapsid (N) gene. This study is the first report on the genetic 
relationship between one systemic (FRSCV MSU-1) and one enteric 
(FRECV MSU-2) ferret coronavirus strain, based upon a comparative 
analysis of their 3’-most 8.3-8.6 kb genomic sequences. 


2. Materials and methods 
2.1. Samples and strain designations 


Tissue samples were obtained from a ferret that had exhibited 
clinical signs and lesions of a FIP-like disease. Immunohistochemi- 
cal staining of the affected tissues (lymph nodes, lung, liver, spleen, 
and brain) with anti-feline coronavirus antibody was positive. Total 
RNA was extracted from the tissue specimens using the QIAGEN 
RNeasy Mini Kit (Valencia, CA). The systemic ferret coronavirus ana- 
lyzed in this study is hereafter referred to as strain FRSCV MSU-1. 

Fecal samples from a ferret diagnosed with ECE were obtained. 
RNA was extracted and was found positive for FRECV by an N-gene 
specific RT-PCR (Wise et al., 2006). The enteric ferret coronavirus 
strain was designated FRECV MSU-2. 


2.2. Deriving the sequences of the distal one-third of the genomes 
of FRSCV MSU-1 and FRECV MSU-2 


The primer sequences pertaining to these experiments are pre- 
sented in Supplementary Material Table 1. Initially, partial S, M 
and N sequences were obtained for each ferret coronavirus strain 
using coronavirus consensus RT-PCR primers (S primers 55 and 
56, and M-N primers 25 and 24 (Tobler and Ackermann, 1996), 
as previously described (Wise et al., 2006) The partial S sequence 
(approximately 600 bases) is located less than 200 nucleotides from 
the 3’ end of the gene, while the M-N target consists of approx- 
imately the last 350 nucleotides of the M gene extending into 
approximately the first 350 nucleotides of the N gene. 

To bridge the approximately 2 kb gap between the 3’ end of the 
S gene and the 3’ end of the M gene, primers were designed from 
partial sequences of the S and M genes. The QIAGEN OneStep RT- 
PCR Kit was used with a0.5 pM final primer concentration ina50 pl 
reaction volume. Cycling conditions were as follows: RT step at 48 C 
for 30 min, hotstart denaturation at 95 C for 15 min, followed by 40 
cycles of 94C for 30s, 55 C for 30s and 72 for 2 min. Final extension 
was at 72C for 7 min. 

The rest of the N gene plus downstream sequence up to the 
3’ untranslated region (UTR) of the genome was obtained by per- 
forming a 3’ RACE, as previously described (Wise et al., 2006). 
Gene-specific forward primers FRECV N-GSP and FRSCV N-GSP 
were designed from available partial 5’ end sequence of the N genes. 

For FRSCV MSU-1, an additional 3’ RACE product that starts from 
the 3’ end of the S gene, was obtained using FRSCV S-GSP as the for- 
ward primer. Complementary DNA was synthesized according to 
the manufacturer’s protocol (3’ RACE System, Invitrogen, Carlsbad, 
CA). PCR of the target cDNA was performed using the Accuprime 
Taq DNA Polymerase High Fidelity Kit (Invitrogen). The 50 jl reac- 
tion mix contained 2 jl of the cDNA template, 0.2 wM of S-GSP and 
AUAP primers, 5 wl of 10x Buffer II, 1 wl of Taq and proofreading 
enzymes, and sterile water. Cycling conditions were as follows: hot- 


start denaturation at 94C for 2 min, followed by 35 cycles of 94C 
for 15s,55C for 30s, and 68 C for 6 min. Final extension was at 68 C 
for 7 min. 

To derive the remaining S gene sequence (~3.2 kb) upstream 
of the already sequenced 3’ end of the gene, 3 overlapping PCR 
fragments (S Codehop 1, S Codehop 2 and Pol 1b-S Codehop 3) were 
amplified sequentially by designing degenerate forward primers 
via the CODEHOP (consensus-degenerate hybrid oligonucleotide 
primer) strategy (Rose et al., 2003) and specific reverse primers 
from already determined sequence data. The CODEHOP degenerate 
PCR primers were obtained by identification of conserved blocks 
of amino acids between the S and pol1b protein sequences of the 
group 1 animal coronaviruses FIPV, CCV and TGEV. The QIAGEN 
OneStep RT-PCR was used for amplifying each of the 3 fragments, 
using 0.5 wM final primer concentrations in 50 wl reaction volumes. 
It was necessary to add the Q solution for S Codehop 2 amplification. 
RT and pre-denaturation steps of 50 C for 30 min and 95 C for 15 min 
were carried out for all 3 fragments. This was followed by 10 cycles 
of touchdown PCR with gradually decreased annealing from 60 C to 
55 C(for S Codehop 1), from 63 C to 56 C (for S Codehop 2), and from 
60 C to 51 C (for Replicase-S Codehop 3), followed by an additional 
of 35 cycles of 94C for 30s, the respective annealing temperatures 
after touchdown for 30s and extension at 72 C for 2-3 min. Final 
extension was at 72 C for 7 min. 

The PCR products were extracted from agarose gels using the 
QIAquick Gel Extraction Kit (QIAGEN). Purified PCR products were 
submitted to the Research Technology Support Facility of Michi- 
gan State University for automated sequencing with both forward 
and reverse amplification primers. This was followed by primer- 
walking with newly designed 5’ and 3’ primers derived from 
acquired sequence data. All specific primers were designed using 
the OLIGO 6 primer analysis software (Molecular Biology Insights, 
Cascade, CO). 

Verification of the first set of sequence data was done by 
direct sequencing of 6 overlapping fragments (~1.5-2 kb each) 
using newly designed sequence-specific primers (Supplementary 
Material Table 2). The final sequence data have an average coverage 
of 4x. 


2.3. Sequencing of the partial spike genes of additional FRSCV and 
FRECV strains 


Samples from additional ECE and FIP-like case submissions were 
used to identify and analyze additional strains. Consensus RT-PCR 
for the S gene was used to derive the partial gene sequences 
of 2 additional FRSCV strains (FRSCV MSU-S and FRSCV WADL) 
and 1 FRECV strain (FRECV 1202). FRSCV WADL was the same 
strain from which the S gene was partially sequenced in a pre- 
vious study (Garner et al., 2008). The partial S sequences obtained 
from these additional strains were compared to the correspond- 
ing sequences of FRSCV MSU-1 and FRECV MSU-2 (this study) 
and FRECV MSU-1 (previously designated as FECV MSU1; GenBank 
Accession DQ340561; Wise et al., 2006). 


2.4. Sequencing of the ORF 3 regions of additional FRSCV and 
FRECV strains 


ORF 3 region sequences were also obtained from strains FRECV 
1202, FRSCV MSU-S and FRSCV WADL, and compared to the 
sequences of FRECV MSU-2 and FRSCV MSU-1. The ORF 3 region 
was amplified using the S-E (spike-envelope) region primer set, 
targeting a 0.8-1kb region (Supplementary Material Table 1). 
Amplification was performed following the same protocol used for 
the S-M region RT-PCR, except for an annealing temperature of 52 C 
and an extension time of 1 min and 15s. 
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2.5. Sequence and phylogenetic analyses 


BLAST (Altschul et al., 1990) similarity searches against the 
GenBank database were performed to determine the nature of 
the sequence obtained. Sequence assembly and analyses, includ- 
ing open reading frame (ORF) searches and multiple alignments of 
nucleotide and predicted amino acid sequences using the Clustal W 
(Thompson et al., 1994), were done with the Lasergene biocomput- 
ing software (DNASTAR, Inc., Madison, WI). Phylogenetic trees were 
constructed with the Treecon software package (Van de Peer and De 
Wachter, 1994), using the neighbor-joining (Saitou and Nei, 1987) 
and 100 bootstrap analyses (Efron and Gong, 1983; Felsenstein, 
1985). Sequence alignment data generated with Lasergene were 
converted to the Treecon format using the ForCon software (Raes 
and Van de Peer, 1999). 


2.6. Sequences and GenBank accession numbers 


The following are the GenBank accession numbers of the 
sequences used in the construction of the phylogenetic trees: 
DQ010921 (FIP 79-1146); EU186072 (FCoV Black); D13096 (CCV 
INSAVC); DQ112226 (CCV CB/05); DQ340562 (FRECV MSU-1); 
EF192159, EF192155 and EF192157 (RDCoV pol 1b, S andN, respec- 
tively); EF192160, EF192156, and EF192158 (CFBCoV pol 1b, S and 
N, respectively); DQ811786 (TGEV); DQ811787 (PRCV); NC_003436 
(PEDV); NC_002645 (HCoV 229E); EU371564 (SARS); GU338456 
(FRSCV MSU-1, this study); GU338457 (FRECV MSU-2, this study). 

The partial S and complete ORF 3 gene sequences obtained from 
additional FRSCV and FRECV strains were deposited in GenBank 
under accession numbers GU459055-GU459060. 


2.7. Development of two genotype-specific RT-PCR assays for the 
detection of ferret coronaviruses in clinical samples 


A multiple alignment of the partial S nucleotide sequences 
(approximately 600 bases) of FRSCV and FRECV strains, correspond- 


ing to the amino acid alignment in Fig. 1, was performed (data 
not shown). The sequences were highly conserved among strains 
within a pathotype group. The 3 FRSCV strains shared an average 
similarity of 94.7%, while the 3 FRECV strains had a 91.8% average 
identity. Average similarity between enteric and systemic strains 
was only 78.5%. These observations made it possible to design two 
genotype-specific RT-PCR primer sets. Fig. 2 shows the aligned 
PCR target regions (sense strand) within the S genes of FRSCV 
MSU-1 and FRECV MSU-2. The positions of the genotype-specific 
primers for each virus are underlined. The genotype 1-specific 
primers, forward, 5’-CTGGTGTTTGTGCAACATCTAC-3’, and reverse, 
5’-TCTATTTGCACAAAATCAGACA-3’, are based upon the FRSCV 
sequences, and amplify a 157-bp product. The genotype 2-specific 
primers, forward, 5’-GGCATTTGTTTTGATAACGTTIG-3’, and reverse, 
5’-CTATTAATTCGCACGAAATCTGC-3’, are based upon the FRECV 
sequences, and amplify a 146-bp product. These primer sets were 
used at a final concentration of 0.5 uM in a 50 pl reaction volume, 
using the QIAGEN One-Step RT-PCR Kit. The RT-PCR assays were 
run with the following cycling parameters: an RT step of 50C for 
30 min, pre-denaturation at 95C for 15 min; followed by 40 cycles 
of 94C for 30s, 53C for 30s and 72C for 30s. Final extension was 
at 72C for 7min. The products were analyzed in 2% pre-stained 
agarose gels and visualized by UV transillumination. 


3. Results 


3.1. Characterization of the sequences spanning the distal 
one-third of the genomes of FRSCV MSU-1 and FRECV MSU-2 


Sequences of the distal third of the genomes of systemic ferret 
coronavirus, FRSCV MSU-1 (8.3 kb), and enteric ferret coronavirus, 
FRECV MSU-2 (8.6kb), were derived. Each sequence spans the 
last 317 nucleotides of the polymerase 1b gene, the full coding 
sequences of structural protein genes (S, E, M and N) and the non- 
structural protein (nsp) genes 3b or 3c-like (ORF 3), 3x-like and 7b 
(ORF7). Partial sequence was also derived from the 3’ untranslated 
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Fig. 1. Alignment of the deduced partial amino acid sequences of the spike proteins of FRSCV MSU-1, FRSCV MSU-S, FRSCV WADL, FRECV MSU-2, FRECV-MSU1 and FRECV 


1202. Residues that match the majority sequence are shown as dots. 
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FRSCV MSU-1 
FRECV MSU-2 eer’ 
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Fig. 2. Alignment of the genotype-specific diagnostic PCR target regions within the spike genes of FRSCV MSU-1 and FRECV MSU-2. Matching nucleotides are shown as dots. 


The broken line indicates a deletion. Primer sequences are underlined. 


region (UTR) of the genome. The genomic organization observed for 
each virus was consistent with that of a typical coronavirus. Each 
intact gene was preceded, a few nucleotides upstream its 5’ end, 
by the recognized consensus transcription-regulating sequence 
(TRS), 5’-CTAAAC-3’, of coronaviruses (Budzilowicz et al., 1985; 
Horsburgh et al., 1992; Kapke and Brian, 1986; Snijder et al., 2003; 
Spaan et al., 1988). In both viruses, the TRS for the E and 7b genes 
each exhibited a one-base mismatch at the 5’ end, appearing as 
5’-GTAAAC-3’ for the former, and 5’-TTAAAC-3’ for the latter. The 
highly conserved 10-base sequence, 5’-GGGAAGAGCT-3’, found in 
the 3’ UTR of all coronaviruses, was also noted, both in FRSCV 
and FRECV. Reminiscent of the observation made for FRECV MSU- 
1 (Wise et al., 2006), the 10-base sequence was found 170 bases 
downstream ORF 7b and had the same 2 base mismatches, a T 
instead of a Gin the first base, and a G instead of an A at the seventh 
base. 

Each ORF was analyzed by BLAST against the GenBank database 
and was confirmed to be of coronavirus origin. Table 1 shows the 
summary of the ORFs identified for each virus strain, the gene and 
predicted protein sequence lengths, and the nucleotide and amino 
acid sequence similarities between FRSCV MSU-1 and FRECV MSU- 
2. Among the nsp genes identified, a single open reading frame, ORF 
3, between the S and E genes was found for both viruses. BLAST 
analysis of the putative ORF 3 protein showed highest significant 
similarities (56-59%) to nsp 3c of CCV and FCoV, and to the 3c homo- 
logue, nsp 3b of TGEV. Two more putative nsp genes observed were 
the two coronavirus ORF 7-like genes downstream of the N gene, 
namely, 3x-like and 7b, as previously observed in the FRECV MSU-1 
Strain (Wise et al., 2006). However, FRSCV MSU-1 only contained 
a 45-bp remnant of the 3x-like ORF because of a 113-bp deletion, 
33 bases downstream of the start codon, resulting to a frameshift 
that prematurely ends the ORF. FRECV MSU-2 shared very high 
sequence similarities to previously reported partial S and M, full 
N, 3x-like and 7b protein sequences of FRECV MSU-1 (96.4-100%). 
FRSCV MSU-1 shared 94.6-100% similarities with FRECV MSU-1 
for the M, N and 7b sequences, but exhibited a marked sequence 
difference (only 77.4% identity) for the partial S protein sequence. 

The ORF lengths were identical for the pol, E, M, N and 7b genes 
of FRSCV MSU-1 and FRECV MSU-2. The 4.35kb S sequence of 
FRECV MSU-2 is 24 nucleotides shorter, observed as a deletion of 
8 amino acid residues in a pairwise comparison of predicted pro- 
tein sequences (Fig. 3). ORF 3 of FRSCV MSU-1 is only 252 bases 
(83 amino acids) long compared to 744 bases (247 amino acids) 
in FRECV MSU-2. In the ORF 3 sequence alignment in Fig. 4, it can 
be observed that the ATG start codon at position 9, found in the 
enteric strain, has been replaced with a GTG in the systemic strain, 
resulting in a delayed start site further downstream at nucleotide 
position 311 for FRSCV MSU-1. In addition, the systemic strain has 
a 193-base deletion that corresponds to positions 228-420 within 
the intact ORF 3 gene of FRECV MSU-2. 

Pairwise alignment of the derived partial genomic sequences of 
FRSCV MSU-1 (8.3 kb) and FRECV MSU-2 (8.6kb) showed only a 
similarity of 87%. For the individual ORFs, the genes that exhibited 
the highest nucleotide sequence similarities (96.2-98.1%) between 
the two viruses were the partial pol1b, full M, N and 7b genes. Sim- 


ilarities based on the deduced amino acid sequences for these 4 
genes ranged from 94.6 to 98.1%. Similarity of the truncated puta- 
tive ORF 3 nsp of FRSCV MSU-1 to the corresponding sequence 
in FRECV MSU-2 was 96.3%. The E gene exhibited intermediate 
nucleotide (91.6%) and amino acid (87.8%) similarities. However, 
significant divergence was observed in the S gene, with nucleotide 
and amino acid sequence similarities of only 79.5 and 79.6%, respec- 
tively. The alignment of the deduced S protein sequences of FRSCV 
MSU-1 and FRECV MSU-2 (Fig. 3) shows the amino acid sub- 
stitutions to be randomly distributed throughout with up to 12 
consecutive substitutions noted in the middle third of the sequence. 


3.2. Analysis of partial spike gene sequences of additional 
systemic and enteric ferret coronaviruses 


To determine whether the marked difference (79.6% amino acid 
similarity) observed in the S genes of FRSCV MSU-1 and FRECV 
MSU-2 may possibly be linked to pathotype difference, and not to 
random heterogeneity of the S sequence in ferret coronaviruses, 
a limited number of additional strains from different geographic 
locations were analyzed. The partial S sequences of two other sys- 
temic strains, FRSCV MSU-S and FRSCV WADL, and one other enteric 
strain, FRECV 1202, were determined by consensus S gene RT-PCR. 
Their deduced amino acid sequences were aligned with corre- 
sponding sequences of FRSCV MSU-1, FRECV MSU-2 and FRECV 
MSU-1 (Fig. 1). The partial sequence being analyzed corresponds 
to residues 1195-1393 of the FRSCV MSU-1 S protein in Fig. 3. The 
alignment in Fig. 1 shows that 21 amino acid substitutions observed 
between the systemic and enteric viruses were conserved within 
a virus pathotype. Moreover, the same four residue deletion was 
consistently observed in each of the three FRECV strains exam- 
ined in comparison to the three strains of FRSCV. The systemic 
strains shared an average identity of 94.3% with each other. The 
average similarity shared between two enteric strains was 92.3%. 
Average similarity of the partial spike sequences between the two 
pathotypes was only 80.1%. 


3.3. Phylogenetic analyses 


The deduced amino acid sequences of the partial polymerase, 
full S, E, M and N proteins of FRSCV MSU-1 and FRECV MSU-2 
were aligned with corresponding sequences of previously reported 
eroup 1 coronaviruses. The phylogenetic trees derived are shown 
in Fig. 5. Sequence similarities between ferret coronavirus and rep- 
resentative group 1 coronaviruses ranged from 66.3 to 74% for the 
partial pol, 43 to 63.6% for S, 32.5 to 64.2% for E, 47.3 to 66.7% for 
M, and 34.6-60.2% for N. The two ferret viruses were observed to 
segregate into their own distinct cluster within the group. 


3.4. Analysis of ORF 3 region sequences of systemic and enteric 
ferret coronaviruses 


ORF 3 region sequences of additional strains, FRECV 1202, FRSCV 
WADL, and FRSCV MSU-S, were determined and compared to the 
sequences of FRECV MSU-2 and FRSCV MSU-1. Enteric strain FRECV 
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Table 1 
Comparison of FRSCV MSU-1 and FRECV MSU-2 sequences spanning the distal one-third 


Gene/protein name ORF length in bp No of amino acids 
FRSCV FRECV FRSCV FRECV 

Pol 1b? SA 7 B47 104 104 
s 4374 4350 1457 1449 
ORF 3 D5 744 83 247 
E 249 249 82 82 
M 792 792 263 263 
N 1125 1125 374 374 
3x-like 45 225 14 74 
7b 615 615 204 204 


nd=not done. 
4 Partial sequence only. 
b FRSCV sequence compared to corresponding partial sequence of FRECV. 


1202 and systemic strain FRSCV MSU-S both have an intact ORF 3, 
consisting of 744 nucleotides, similar to that of FRECV MSU-2. Sim- 
ilarity of the intact putative ORF 3 protein between the two enteric 
strains was only 79.8%. The ORF 3 protein of FRSCV MSU-S was only 
78.1% similar to FRECV 1202, but shared a much higher identity 
with enteric strain, FRECV MSU-2, at 93.9%. The ORF 3 sequence of 


rch 149 (2010) 42-50 


of the genome. 


% Amino acid similarity between 
FRSCV and FRECV 


% Nucleotide similarity between 
FRSCV and FRECV 


98.1 98.1 
79.6 79.5 
96.3> 96.4» 
87.8 91.6 
97 97.7 
94.9 96.2 
nd nd 
94.6 96.6 


systemic strain FRSCV WADL had a deletion of 86 bases (positions 
256-341 of FRECV MSU-2 in Fig. 4) that overlaps the 193-base dele- 
tion in FRSCV MSU-1. This deletion results to the splitting of ORF 3 
into two shorter ORFs in FRSCV WADL—a proximal ORF, express- 
ing a putative protein of 98 amino acids with a start site similar 
to that of FRECV MSU-2, and a distal ORF that translates into a 
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Fig. 3. Alignment of the deduced amino acid sequences of the spike proteins of FRSCV MSU-1 and FRECV MSU-2. Conserved residues are shown as dots. 
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Fig. 4. Alignment of the ORF 3 region sequences of FRECV MSU-2 and FRSCV MSU-1. Matching nucleotides are shown as dots. Deletions are represented by broken lines. ATG 


and TAG or TAA codons are underlined. 


110 amino acid protein, whose ATG corresponds to position 420 in 
FRECV MSU-2 (Fig. 4). ORF 3 region sequences of systemic strains 
FRSCV WADL and MSU-S were highly similar at 95.5%. The corre- 
sponding ORF 3 region sequence of FRSCV MSU-1 appeared more 
distant with an average identity of only 85.75% to either systemic 
Strain. These results indicate high sequence variability in the ORF 3 
region of the ferret coronavirus strains examined. 


3.5. Detection of FRSCV or FRECV in specimens using the S-gene 
genotype-specific diagnostic RT-PCR primer sets 


RNA from various tissue samples, from which strains FRSCV 
MSU-1, FRSCV-WADL and FRSCV MSU S had been previously iden- 
tified, were tested with the genotype 1-specific diagnostic primer 
set. The 157-bp S gene target was amplified from liver, lung, spleen, 
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Fig. 5. Phylogenetic trees based on deduced partial amino acid sequences of the: (A) partial polymerase, (B) spike, (C) envelope, (D) membrane, and (E) nucleocapsid proteins 
of FRSCV MSU-1, FRECV MSU-2 and representative strains of group 1 coronaviruses. SARS virus was included as the outgroup sequence. Bootstrap values are indicated at the 


nodes. 
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Fig. 6. S-gene genotype-specific diagnostic RT-PCR on clinical samples. (A) Geno- 
type 1-specific assay, 157-bp positive amplicon observed only on Lanes 2-6. (B) 
Genotype 2-specific assay, 147-bp positive amplicon observed only on Lanes 7-8. 
Lane designations for both panels A and B: Lanes 1 and 10, 100-bp DNA ladder; Lane 
2, FRSCV MSU-1 positive lung; Lane 3, FRSCV MSU-1 positive kidney; Lane 4, FRSCV 
WADL positive lymph node; Lane 5, FRSCV MSU-S positive spleen; Lane 6, FRSCV 
MSU-S positive intestine; Lane 7, FRECV MSU-2 positive feces; Lane 8, FRECV 1202 
positive feces; Lane 9, negative control (sterile water). 


kidney, brain, lymph node and intestinal tissues collected from fer- 
rets affected with the FIP-like disease. Fig. 6A shows the specific 
157-bp product amplified from each sample for a subset of these 
specimens. Enteric strains FRECV MSU-2 and FRECV 1202, did not 
amplify with the genotype 1-specific primers. Fig. 6B shows the 
results for the same set of samples tested in Fig. 6A, but using 
the genotype 2-specific primer set. The genotype 2-specific RT-PCR 
assay detected the 147-bp S gene target only in the samples known 
to contain FRECV strains. 


4. Discussion 


Ferret systemic coronavirus causes signs and lesions that are 
similar to those induced by FIPV in cats. The partial S protein 
sequence we initially derived from a FRSCV strain (Garner et al., 
2008), was the first indicator that this virus appeared to be more 
similar to ferret enteric coronavirus than to any other group 1 
coronavirus. The question then arose as to whether FRSCV strains 
could have evolved from FRECV by in vivo mutation, leading to the 
acquisition of a systemic virulent pathotype. The analogous situa- 
tion in cats is enteric feline coronavirus infection, which causes a 
mild or asymptomatic localized enteric infection (Pedersen, 2009; 
Evermann et al., 1991; Rottier, 1999; Weiss and Navas-Martin, 
2005; Addie et al., 2004) and the systemic highly fatal disease FIP, 
considered to be derived by spontaneous mutation during FCoV 
persistence (Evermann etal., 1991; Rottier, 1999; Weiss and Navas- 
Martin, 2005; Addie et al., 2004; De Groot-Mijnes et al., 2005; 
Poland et al., 1996; Vennema et al., 1998; Herrewegh et al., 1995; 
Rottier et al., 2005; Lai et al., 2007; Pedersen, 2009; Chang et al., 
2009). The mutation from FCoV to FIPV is accompanied by a marked 
acquisition of macrophage tropism (Pedersen, 1987; Stoddart and 
Scott, 1989), which facilitates systemic spread (Rottier et al., 2005; 
Lai et al., 2007; Pedersen, 1987). The precise nature of this phe- 
nomenon is not entirely understood (Lai et al., 2007). In fact, enteric 
FCoVs and their counterpart FIPVs were observed to be geneti- 


cally very similar (Poland et al., 1996; Herrewegh et al., 1995; 
Vennema et al., 1998), bearing only minor but crucial differences 
(Pedersen, 2009). Candidate genes speculated to play a role in the 
virulence shift include the S and group-specific genes 3abc, 7b 
and 7a (Haijema et al., 2004; Vennema et al., 1998; Kennedy et 
al., 2001). Using targeted RNA recombination, Rottier et al. (2005) 
have shown that the C-terminal 582-amino acid region of the FIPV 
79-1146 S protein contains the amino acid changes that determine 
macrophage tropism. In this region of the S protein, a total of only 
10 amino acid differences were noted between FIPV 79 and 1146 
and its close relative, enteric strain FCoV 79-1683. A “D-to-A” sub- 
stitution, upstream of the putative fusion peptide, was speculated 
to play a role in the transition of FCoV to FIPV. Pedersen (2009), 
in a comprehensive review of FIPV infection, states that mutation 
in the 3c gene is responsible for the FIP biotype. Similarly, Chang 
et al. (2009) analyzed the 3c gene in 27 feline enteric coronavirus 
(FECV)- and 28 FIPV-infected cats and found that the 3c gene was 
mutated in the majority (71.4%) of the cats with FIP. 

An alternative to the in vivo mutation hypothesis of FIPV strain 
development was put forward by Brown et al. (2009) and was based 
on phylogenetic analyses of the M and nsp 7b gene sequences from 
FCoV-infected clinically healthy and sick cats. They observed that 
viral sequences from healthy cats were distinct from those of sick 
cats supporting the hypothesis that distinctive virulent and aviru- 
lent FCoV strains co-circulate in natural populations. 

With the currently available data, it is still difficult to predict 
whether the systemic ferret coronavirus, FRSCV, is derived from 
its enteric counterpart, FRECV, by in vivo mutation, or is a co- 
circulating distinct strain. To resolve this, a significant number of 
additional strains from the same or different locations will have to 
be examined. It is also not clear at this point whether FRSCV orig- 
inated as a result of a recombination between FRECV and another 
group | coronavirus, resulting in the acquisition of a very distinct S 
gene. A 3’/-5’ scanning of the nucleotide sequence alignment for the 
two virus strains did not show an abrupt discontinuity in identity 
at the 3’ end of the S gene that would have been an indication of 
recombination. We noticed instead an increase in sequence diver- 
gence Starting at the 200 bases adjacent to the 3’ end of the S gene, 
within the ORF 3 gene region (nucleotide positions 1-200 in Fig. 4), 
that gave a sequence similarity of only 75.5%. This of course does 
not discount the possibility that a recombination event may have 
involved the proximal portion of the ORF 3 gene region as well. 
BLAST analysis of the S sequence of FRSCV MSU-1 against reported 
sequences only showed FRECV MSU-1, the previously described 
enteric ferret coronavirus strain (Wise et al., 2006), as the one with 
the most significant similarity to it. If recombination in the S gene 
did occur, giving rise to FRSCV, the donor coronavirus strain is yet 
to be discovered. An example of this phenomenon is type II FCoV, 
which originated from a double recombination between FCoV type 
I and canine coronavirus (CCV) (Motokawa et al., 1996; Herrewegh 
et al., 1998). 

The main conclusion that can be reached from the present data 
is that there is a significant difference observed between the S genes 
of FRSCV MSU-1 and FRECV MSU-2. Whether this difference is asso- 
ciated with pathotype is not clear. In vivo experiments would be 
an obvious way to resolve this issue, but are hampered by the diffi- 
culty of growing either virus in cell culture. Another observation 
was that the partial S gene sequences from at least 3 different 
FRSCV strains, collected from different locales, were distinct from 
3 different FRECV strains that also originated from different places 
(Fig. 1). Whether this pattern would hold up with a larger sample 
size remains to be determined. The two genotype-specific RT-PCR 
assays developed in this study would certainly be useful molecular 
tools to expand the current data. 

Another proposed virulence marker in FIP strains is the 3c gene 
(Vennema et al., 1998). According to Pedersen (2009), mutation in 
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the ORF 3c gene is responsible for the FIPV biotype and that FCoVs 
of the enteric biotype possess an intact 3c gene. It is thought that a 
loss of 3c gene function enhances the internalization and replica- 
tion of FIPVs in macrophages (Dewerchin et al., 2005). Pedersen 
(2009) explains that the mutations could either be SNPs which 
lead to premature stop codons or deletion mutations that signif- 
icantly truncate the 3c gene product. In agreement to this are the 
findings of Chang et al. (2009) where mutations in the 3c gene 
were found in a majority (71.4%), of FIPV strains examined, while 
the gene was intact in 100% of analyzed FECVs. This lead to the 
conclusion that mutation in 3c cannot be the sole cause of FIP. 
In the ferret strains we analyzed, both systemic and enteric bio- 
types possess a single ORF 3 protein. The 3c-like protein in FRSCV 
MSU-1 was significantly truncated in comparison to FRECV MSU-2, 
prompting us to explore other systemic and enteric strains. Our 
findings show that the ferret coronaviruses have high sequence 
variability in ORF 3, even among strains of the same pathotype. 
A deletion in the 3c-like gene, giving rise to a truncated protein, 
was noted in two out of three systemic strains examined, while an 
intact ORF 3 was noted in the two enteric strains analyzed. Again, 
with the limited number of strains studied so far, it is not clear 
whether the mutation in gene 3 plays a role in the systemic nature 
of FRSCVs. 

In conclusion, this study presents the first molecular compar- 
ison of two ferret coronavirus strains of different pathotypes. We 
provide molecular evidence showing that systemic strain, FRSCV 
MSU-1, is more closely related to enteric strain, FRECV MSU-2, 
than to any other group 1 coronavirus, in all structural and nsp 
genes examined. One important observation was that these two 
virus strains differed significantly in the S gene. It is conceivable 
that FRSCV and FRECV are genetically distinctive coronaviruses of 
ferrets, a hypothesis already supported by Brown et al. (2009) in 
the case of FCoVs. On the other hand, it is also possible that the in 
vivo mutation theory of Chang et al. (2009) and Pedersen (2009) 
applies to the ferret coronaviruses. Analysis of more strains of each 
virus pathotype is certainly needed to either substantiate or modify 
either of these two concepts. 


Acknowledgments 


We thank Ross Burwell for the excellent technical help, and 
the DCPAH Virology staff for assisting with sample processing and 
tracking. We also thank Dr. Jim Evermann of the Washington Ani- 
mal Disease Diagnostic Laboratory at Pullman, WA, for sending 
ferret diagnostic tissue samples. This study was funded by a grant 
from the Companion Animal Fund of the College of Veterinary 
Medicine, Michigan State University. 


Appendix A. Supplementary data 


Supplementary data associated with this article can be found, in 
the online version, at doi:10.1016/j.virusres.2009.12.011. 


References 


Addie, D.D., Paltrinieri, S., Niels, C.P., 2004. Recommendations from workshops of the 
second international feline coronavirus/feline infectious peritonitis symposium. 
J. Feline Med. Surg. 6, 125-130. 

Altschul, S.F., Gish, W., Miller, W., Myers, E.W., Lipman, D.J., 1990. Basic local align- 
ment search tool. J. Mol. Biol. 215, 403-410. 

Brown, M.A., Troyer, J.L., Pecon-Slattery, J., Roelke, M.E., O’Brien, S.J., 2009. Genetics 
and pathogenesis of feline infectious peritonitis virus. Emerg. Infect. Dis. 15, 
1445-1452. 

Budzilowicz, C.J., Wilczynski, S.P., Weiss, S.R., 1985. Three intergenic regions of 
coronavirus mouse hepatitis virus strain A59 genome RNA contain a common 
nucleotide sequence that is homologous to the 3’ end of the viral mRNA leader 
sequence. J. Virol. 53, 834-840. 

Chang, H.W., de Groot, R.J., Egberink, H.F., Rottier, P.J.M., 2009. Feline infectious 
peritonitis: insights into feline coronavirus pathobiogenesis and epidemi- 


ology based on genetic analysis of the viral 3c gene. J. Gen. Virol. 
doi:10.1099/vir.0.016485-0. 

De Groot-Mijnes, J.D.F., Van Dun, J.M., Van der Most, R.G., De Groot, R.J., 2005. Nat- 
ural history of a recurrent feline coronavirus infection and the role of cellular 
immunity in survival and disease. J. Virol. 79, 1036-1044. 

Dewerchin, H.L., Cornelissen, E., Nauwynck, H.J., 2005. Replication of feline coron- 
aviruses in peripheral blood monocytes. Arch. Virol. 150, 2483-2500. 

Efron, B., Gong, G., 1983. A leisurely look at the bootstrap, the jackknife, and cross- 
validation. Am. Stat. 37, 36-48. 

Evermann, J.F., McKeirnan, A.J., Ott, R.L., 1991. Perspectives on the epizootiology of 
feline enteric coronavirus and the pathogenesis of feline infectious peritonitis. 
Vet. Micro. 28, 243-255. 

Felsenstein, J., 1985. Confidence limits on phylogenies: an approach using the boot- 
strap. Evolution 39, 783-791. 

Garner, M.M., Ramsell, K., Morera, N., Juan-Sallés, C., Jiménez, J., Ardiaca, M., Mon- 
tesinos, A., Teifke, J.P., Lohr, CV., Evermann, J.F., Baszler, T.V., Nordhausen, R.W., 
Wise, A.G., Maes, R.K., Kiupel, M., 2008. Clinicopathologic features of a systemic 
coronavirus-associated disease resembling feline infectious peritonitis in the 
domestic ferret (Mustela putorius). Vet. Pathol. 45, 236-246. 

Haijema, B.J., Volders, H., Rottier, P.J.M., 2004. Live, attenuated coronavirus vaccines 
through the directed deletion of group-specific genes provide protection against 
feline infectious peritonitis. J. Virol. 78, 3863-3871. 

Herrewegh, A.A.P.M., Vennema, H., Horzinek, M.C., Rottier, P.J.M., de Groot, R.J., 1995. 
The molecular genetics of feline coronaviruses: comparative sequence analysis 
of the ORF7a/7b transcription unit of different biotypes. Virology 212, 622-631. 

Herrewegh, A.A.P.M., Smeenk, I., Horzinek, M.C., Rottier, P.J.M., De Groot, R.J., 1998. 
Feline coronavirus type II strains 79-1683 and 79-1146 originate from a double 
recombination between feline coronavirus type I and canine coronavirus. J. Virol. 
72, 4508-4514. 

Horsburgh, B.C., Brierley, I., Brown, T.D.K., 1992. Analysis of a 9.6 kb sequence from 
the 3’ end of canine coronavirus genomic RNA. J. Gen. Virol. 73, 2849-2862. 
Juan-Sallés, C., Teifke, N., Morera, N., Jimenez, J., Montesinos, A., Ardiaca, M., Lohr, 
C.V., Garner, M.M., 2006. Pathology and immunohistochemistry of a disease 
resembling feline infectious peritonitis in ferrets (Mustela putorius furo). Proc. 

Am. Col. Vet. Pathol. 84, 845. 

Kapke, P.A., Brian, D.A., 1986. Sequence analysis of the porcine transmissible gas- 
troenteritis coronavirus nucleocapsid protein gene. Virology 151, 41-49. 

Kennedy, M., Boedeker, N., Gibbs, P., Kania, S., 2001. Deletions in the 7a ORF of feline 
coronavirus associated with an epidemic of feline infectious peritonitis. Vet. 
Micro. 81, 227-234. 

Lai, M.M.C., Perlman, S., Anderson, L.J., 2007. Coronaviridae. In: Knipe, D.M., How- 
ley, P.M. (Eds.), Fields Virology, 5th edition. Lippincott Williams & Wilkins, pp. 
1305-1335. 

Martinez, J., Ramis, A.J., Reinacher, M., Perpinan, D., 2006. Detection of feline infec- 
tious peritonitis virus-like antigen in ferrets. Vet. Rec. 158, 523. 

Martinez, J., Reinacher, M., Perpifian, Ramis, A., 2008. Identification of group 1 
coronavirus antigen in multisystemic granulomatous lesions in ferrets (Mustela 
putorius furo). J. Comp. Pathol. 138, 54-58. 

Motokawa, K., Hohdatsu, T., Hashimoto, H., Koyoma, H., 1996. Comparison of the 
amino acid sequence and phylogenetic analysis of the peplomer, integral mem- 
brane and nucleocapsid proteins of feline, canine and porcine coronaviruses. 
Microbiol. Immunol. 40, 425-433. 

Pedersen, N.C., 1987. Virologic and immunologic aspects of feline infectious peri- 
tonitis virus infection. Adv. Exp. Med. Biol. 218, 529-550. 

Pedersen, N.C., 2009. A review of feline enteric coronavirus infection. J. Feline Med. 
Surg. 11, 225-258. 

Poland, A.M., Vennema, H., Foley, J.E., Pedersen, N.C., 1996. Two related strains 
of feline infectious peritonitis virus isolated from immunocompromised cats 
infected with a feline enteric coronavirus. J. Clin. Micro. 34, 3180-3184. 

Raes, J., Van de Peer, Y., 1999. ForCon: a software tool for the con- 
version of sequence alignments. EMBnet.news 6(1). Available at: 
http://www.hgmp.mrc.ac.uk/embnet.news./vol6_1/. 

Rose, T.M., Henikoff, J.G., Henikoff, S., 2003. CODEHOP (Consensus-Degenerate 
Hybrid Oligonucleotide Primer) PCR primer design. Nucl. Acids Res. 31, 
3763-3766. 

Rottier, P.J.M., 1999. The molecular dynamics of feline coronaviruses. Vet. Micro. 69, 
117-125. 

Rottier, P.J.M., Nakamura, K., Schellen, P., Volders, H., Haijema, B.J., 2005. Acquisition 
of macrophage tropism during the pathogenesis of feline infectious peritonitis 
is determined by mutations in the feline coronavirus spike protein. J. Virol. 79, 
14122-14130. 

Saitou, N., Nei, M., 1987. The neighbor-joining method: a new method for recon- 
structing phylogenetic trees. Mol. Biol. Evol. 4, 406-425. 

Snijder, E.J., Bredenbeek, P.J., Dobbe, J.C., Thiel, V., Ziebuhr, J., Poon, L.L.M., Guan, 
Y., Rozanov, M., Spaan, W.J.M., Gorbalenya, A.E., 2003. Unique and conserved 
features of genome and proteome of SARS-coronavirus, an early split-off from 
the coronavirus group 2 lineage. J. Mol. Biol. 331, 991-1004. 

Spaan, W., Cavanagh, D., Horzinek, M.C., 1988. Coronaviruses: structure and genome 
expression. J. Gen. Virol. 69, 2939-2952. 

Stoddart, C.A., Scott, F.W., 1989. Intrinsic resistance of feline peritoneal macrophages 
to coronavirus infection correlates with in vivo virulence. J. Virol. 64, 
436-440. 

Thompson, J.D., Higgins, D.G., Gibson, T.J., 1994. CLUSTAL W: improving the sensi- 
tivity of progressive multiple sequence alignment through sequence weighting, 
position-specific gap penalties and weight matrix choice. Nucl. Acids Res. 22, 
4673-4680. 


50 A.G. Wise et al. / Virus Research 149 (2010) 42-50 


Tobler, K., Ackermann, M., 1996. Identifikation und charakterisierung von neuen 
und unbekannten coronaviren mit hilfe von RT-PCR und degenerierten primern. 
Schweiz. Arch. Tierheilk. 138, 80-86. 

Van de Peer, Y., De Wachter, R., 1994. TREECON for Windows: a software package for 
the construction and drawing of evolutionary trees for the Microsoft Windows 
environment. Comput. Appl. Biosci. 10, 569-570. 

Vennema, H., Poland, A., Foley, J., Pedersen, N.C., 1998. Feline infectious peritonitis 
viruses arise by mutation from endemic feline enteric coronaviruses. Virology 
243, 150-157. 


Weiss, S.R., Navas-Martin, S., 2005. Coronavirus pathogenesis and the emerging 
pathogen severe acute respiratory syndrome coronavirus. Microbiol. Mol. Biol. 
Rev. 69, 635-664. 

Williams, B.H., Kiupel, M., West, K.H., Raymond, J.T., Grant, C.K., Glickman, L.T., 2000. 
Coronavirus-associated epizootic catarrhal enteritis in ferrets. J. Am. Vet. Med. 
Assoc. 217 (4), 526-530. 

Wise, A.G., Kiupel, M., Maes, R.K., 2006. Molecular characterization of a novel coro- 
Navirus associated with epizootic catarrhal enteritis (ECE) in ferrets. Virology 
349, 164-174. 


